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Abstract 

CO 

We discuss a simple approach to solve the transport equation for high-energy neu- 
trinos in media of any thickness. We present illustrative results obtained with some 
specific models for the initial spectra of Vn and Vu propagating through a normal 



*^ cold medium. 
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1 Introduction 



In passing through a medium, neutrinos (and antineutrinos) are absorbed 
and lose their energy due to charged and neutral current interactions. For a 
normal cold medium, like the earth's or stellar interior, these are vN {pN) 
and ve (Ve) collisions. In more exotic media, like the relic black-body neu- 
trino background [1] or hot galactic halos filled by massive neutrinos [2,3], 
ultrahigh energy neutrinos may scatter elastically or be absorbed due to vv 
interactions. Owing to the energy loss and the strong energy dependence of 
total cross sections for neutrino interactions, the neutrino "depth-intensity re- 
lation" (or penetration coefficient [4]) does not follow a simple absorption law 
and the magnitude of this effect grows with energy and depth. Such a situ- 
ation is well known in the muon transport theory (see e.g. ref. [5]). To solve 
the neutrino transport equation for moderate depths, the method of succes- 
sive generations is workable [6]. But this method becomes inefficient for the 
depth in excess of several neutrino interaction lengths. Pertinent refining of the 
neutrino transport theory is desirable for many applications, specifically for 
studying standard and speculative neutrino interactions (see e.g. ref. [7] and 
references therein), detecting neutrino signals from annihilation of dark mat- 
ter particles in the sun and the earth, and for high-energy neutrino astronomy 
with future, km 3 -scale neutrino telescopes [8,9]. 

The goal of this work is to provide an elementary method for the precise cal- 
culation of the energy spectra of high-energy neutrinos after their propagation 
through a medium of any thickness. The problem was considered recently in 
ref. [10] in the framework of a simplified models for the neutrino cross sections 
and initial neutrino spectrum. Our approach does not require simplifications 
and is applicable to cross sections (differential and total) and initial spectra 
of any form. In sections 2 and 3, we shall consider only decreasing unbroken 
initial spectra most interesting for high-energy neutrino astrophysics. How- 
ever, the main idea of the method can be extended also to a monochromatic 
spectrum (see appendix A). This generalization may be of utility, in particu- 
lar, for simulating single neutrino events in a neutrino telescope. Furthermore, 
the method makes no assumptions specific to neutrino transport and can be 
extended almost straightforwardly to the problem of transport of high-energy 
particles other than neutrinos (e.g. hadrons and muons). 

Under certain conditions, neutrinos may transform, changing their flavor via 
processes like V e e~ — > or z/^e~ — > v e l~ (£ ^ e), owing to production 
and decay of short-lived hadrons (D, D s , etc.), or (in the mentioned neutrino 
fields) through reaction chains like v^Vj. — > /i~T + , t + — > T7 T X, etc. As it 
pointed out recently [11], high-energy tau neutrinos (and antineutrinos) will 
effectively regenerate in matter, losing energy, through the charged-current 
reaction chain u T N — > rX, r — > v T X. Such mechanisms must be taken into 
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account in data processing from many future experiments (detecting v T events 
from astrophysical neutrino oscillations at energies > 1 PeV [12], events with 
energy release well beyond the Greisen-Zatsepin-Kuz'min cutoff [3], etc.). 

Mathematically, the inclusion of processes changing the neutrino flavor or 
neutrino energy loss through creation and decay of short-lived particles leads 
to a system of transport equations. The extension of the method to the general 
system is not straightforward and demands additional assumptions specific 
to the task. However, the case when these contributions may be treated as 
corrections presents no special problem. Since this case is rather common (the 
neutrino production in the ^-induced hadronic cascades is a typical example), 
we brief the corresponding trivial generalization in appendix B. 

To avoid technical complications, in the main text we shall neglect the (stan- 
dard and hypothetical) flavor-changing neutrino interactions^ and use the 
simplest "standard" scenario for neutrino propagation described by a single 
transport equation. We shall consider sufficiently high energies in order to ne- 
glect the thermal velocities of the scatterers in the target medium and to deal 
with the one- dimensional theory. As an illustration, we shall discuss results 
obtained with some specific models for the initial spectra of muon neutrinos 
and antineutrinos propagating through a normal cold medium. 



2 Method for the solution of the neutrino transport equation 

Let F u (E,x) be the differential energy spectrum of neutrinos at a column 
depth x in the medium defined by 



where p(L) is the density of the medium at a distance L from the boundary 
measured along the neutrino beam path. Then, neglecting the flavor-changing 
processes mentioned in the introduction, one can derive the following one- 
dimensional transport equation 
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As well as the effects of possible neutrino flavor mixing. 
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with the boundary condition F u (E,0) = F®(E). Here \ V {E) is the neutrino 
interaction length defined by the equation 

jwrV Ta ' S[E) ' 



where N T is the number of scatterers T in 1 g of the medium, affi(E) is the 
total cross section for the vT interactions and the sum is over all scatterer 
types (T = N, e, . . .). The "regeneration function" & u (y, E) is defined by 



where da V T^ v x{y,E)/ dy is the differential cross section for the inclusive re- 
action vT — > vX (with E the initial neutrino energy and y the fraction of 
energy lost) and E y = Ej (1 — y). 

Let us define the effective absorption length A u (E,x) by 

F v (E,x)=F2(E)exp 



x 



K(E,x) 



(2) 



As is evident from eq. (1), A U (E, x) > X V (E) for any finite E and x. Therefore 



where Z v (E,x) is a positive function (we will call it Z factor in analogy with 
the hadronic cascade theory) which contains the complete information about 
the neutrino kinetics in matter. 

Substituting eqs. (2) and (3) into eq. (1) and integrating by parts, it is easy 
to derive the integral equation for the Z factor: 

1 x 1 

Z U (E, x) = -f J V u(y, E)<$> u (y, E) exp [—x'D u (E, E y , x')} dx' dy, (4) 





with 

n ( p p ^ l-Z u (E y ,x) l-Z v (E,x) 

Dv{E ^ x) = ^m xm~ 
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and 



Vu(y,E) = 



FS(E)(i- y y 



We dwell on eq. (4). Although nonlinear, it proves to be more amenable for an 
iteration solution than eq. (1), considering the smallness of the Z factor and 
(what is more important) a model-independent feature of the regeneration 
function § u (y,E), namely, its sharp maximum at y — 0. 

In this section, we assume that the initial spectrum F°(E) is a continuous 
function decreasing at high energies fast enough so that < r] u (y, E) < oo for 
any E and < y < 1. Actually, the neutrino spectra of interest for high-energy 
neutrino astrophysics decrease everywhere so fast that < r) u (y,E) < 1 for 
any E and y > 0. 

We will first look at the case of thin absorbers. One can readily see that 

i 

Z U (E, 0) = / Vu (y, E)<S> u (y, E) dy = Z° U (E). 



The approximation Z v (E,x) = Z®(E) is usually utilized when studying the 
muon neutrino propagation through matter (see e.g. ref. [4,8] and references 
therein). However, at high energies this approximation becomes too rough 
even for "shallow" (as compared to A„) depths. Indeed, taking into account 
the 0(x/X v ) correction yields 



Z„(E,x) 



xAj(E) 
2X U (E) ' 



where 



Al(E) = -X u (E) 



dZ u (E,x) 
dx 



J x=0 



= J v,(y, E)*u(v, e) I [l - Zl{E y )\ ^| - 

Thus, the approximation Z v fa Z® can only be valid for 

2Zl{E) 
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In the general case, the function A*(J5) is not small in comparison with Z®(E). 
This can be demonstrated with the simple model adopted in ref. [10]. The au- 
thors of ref. [10] assumed that & u = &„(y) is an energy independent function, 
X„(E) oc E- p and F°(E) oc E~^ +v> with energy independent positive (3 and 
7. All these assumptions are far from reality but may have a physical sense at 
super-high energies. For example, owing to the vN interactions, /3 is a mono- 
tonically decreasing function of E changing from about 1 at E < 1 TeV to 
about 0.4 at E > 1 PeV [13]; the function ® v (y, E) strongly varies with E at 
all energies, but for E > 1 PeV it may be roughly approximated by a scaling 
function (see fig. 1 in sect. 3).0 In this model, both and A J, are energy 
independent: 



1 

Z° u = J(l- y y<l> u (y)dy, 



Al= (l - Z°) /(l - J/)*" [(1 - y)"" - 1] My) dy. 



Evidently A* can be much larger than Z® for a sufficiently hard initial neutrino 
spectrum (small 7)Q. 



It is not a hard task to derive the O {^x/\ v ) k ^j corrections for k = 2,3, . . ., 
but as a result we will get an asymptotic expansion with coefficient functions, 
A*(E), increasing fast with k. The region of applicability of this expansion 
proves to be very limited and decreases fast with increasing energy. 



Now, let us consider a way to solve eq. (4) for any depth and energy. We will 
use an iteration algorithm. Let n label the iteration number. Then we define 



D<?\E,E v ,x) 



\-Zf\E, 



\ V {Ey 



Zj n) (E,x) 
K{E) 



(5) 



and 



X 1 



Zj? +1 \E, 



x) = — 

X 



Vu {y, E)$ v (y, E) exp \—x'D^ l \E, E y , x')] dx' dy. (6) 







Due to the mentioned sharp maximum of § u {y,E), the main contribution 
into the integral over y on the right side of eq. (6) comes from the lover limit 

2 However, as fig. l.b suggests, the specific parametrization used in ref. [10], <& v {y) oc 
(const + y)^ 1 , is too rough even at super-high energies. 

3 In the real case, A^{E) is nevertheless finite because any physical spectrum F„(E) 
has a cutoff at some finite energy E cnt and therefore y u (y, E) = at y > 1 — E/E cut . 
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neighborhood. But D U (E, E y ,x) — > as y — > and thus the algorithm is 
robust in respect to choosing the zero approximation. The simplest choice is 
Z®(E,x) = 0. Therefore 



The algorithm (5-7) is formally applicable for arbitrary decreasing initial spec- 
tra. It is however clear that the softer the initial spectrum, the better the 
convergence of the algorithm. In the next section, we show that the algorithm 
converges very fast for realistic initial spectra and has no restrictions in depth 
or energy. Moreover, even the first approximation, 



Zi 1 \E,x)= Jr] u (y,E)<S> u (y,E) 



I _ e -xVv{E,E y ) 

xV u {E,E y ) 



(8) 



proves to be quite accurate. It has the correct asymptotic behavior both in en- 
ergy and depth and can thus be used for an analytical or numerical evaluation 
of the Z factor with a not-too-big error. 

Assuming that \ U (E) is a decreasing function^ and therefore V u (E,E y ) > 
for y > 0, one can prove that 

• Zjp(E,x) < Z°(E) at any E and x > 0, and 

• Z^(E, x) — > as x — > oo at any E. 

The latter signifies that neutrino "regeneration" due to the inclusive reactions 
vT — > vX becomes negligible at sufficiently large depths. This conclusion 
remains true for the exact solution of eq. (4), under rather general assumptions 
about the behavior of the initial neutrino spectrum and cross sections. 



3 Numerical Illustration and Discussion 



Below, we will consider only muon neutrinos and antineutrinos propagating 
through a normal medium. But, with obvious reservations, the results that 
follow also hold for electron neutrinos. In this calculation, we shall neglect neu- 
trino scattering off electrons (thus we exclude electron antineutrinos from our 
consideration) as well as the neutrino "recreation" in the //-induced cascades. 
For further simplification, we shall deal with an isoscalar medium neglecting 
nuclear effects. 

4 For a normal medium, this is true for all neutrino flavors except V e (see ref. [13]). 
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To calculate the differential v^N and V^N cross sections we use the approach 
of ref. [13] based on the renormalization-group-improved parton model and 
new experimental information about the quark structure of the nucleon (see 
appendix C). Various versions of different sets of parton density functions 
q(x, Q 2 ) are now collected in a large CERN program library PDFLIB [14]; they 
can be simply accessed by setting few parameters to choose the desired version. 
In this calculations, we selected, following ref. [13], the third version of the 
CTEQ collaboration model [15], because it is characterized by a particularly 
suitable extrapolation at very low Bjorken x. The evolution in Q 2 is realized by 
next-to-leading order Altarelli-Parisi equations from initial Qq = 2.56 GeV 2 . 

The total cross sections for CC and NC inelastic scattering of muon neutrinos 
and antineutrinos off an isoscalar nucleon are shown in fig. l.a as the solid (z/ M ) 
and dashed curves. Fig. l.b shows the regeneration functions $ v (y,E) 
(solid curves) and ^ v (y, E) (dashed curves) versus y for several values of E 
(10 3 to 10 12 GeV). 




Fig. 1. Total (CC and NC) u^N and V^N cross sections vs energy (a) and regen- 
eration functions ^ Ufl (y,E) and <& Vli (y,E) vs y for E = 10 k GeV [k = 3,4, . . . , 12 
from top to bottom] (b). 

At all energies, our calculation for the cross sections agrees with the result of 
ref. [13] within a few percent accuracy; the insignificant difference near the 
resonance region is due mainly to differences in the adopted values for the 
electroweak parameters (W/Z boson masses, t quark mass, Weinberg angle, 
etc.)[[] and, at superhigh energies, to the top sea contribution neglected in 
ref. [13]. As one can see from the figures, the and scatterings become 
indistinguishable for E > 1 PeV. 



5 In our calculation, all these parameters were updated according to the PDG 
data [16]. 
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We use the following model for the initial neutrino spectrum: 



w =fl-('tr7i + £yv#-). 



E J EqJ \ E cnt 



where K, 7, a, E and E cut are parameters and <p(t) is a function equal to at 
t > 1 and 1 at t <C 1. Varying the parameters in eq. (9), we can approximate 
many models for the neutrino fluxes expected from the known astrophysical 
sources. Technically, the function 0(t) serves to avoid an extrapolation of the 
cross sections to the ultrahigh energy region for which our knowledge of the 
parton density functions becomes doubtful. For realistic values of the parame- 
ters 7, a and E , the explicit form of <p(t) is of no importance if one is interested 
in the energy range E <C E cut . In fact, (f>(t) may be treated as the real physical 
cutoff of the spectrum determined by the energetics of the neutrino source or 
by neutrino interactions with the cosmic backgrounds. In the present calcula- 
tions, we adopt (without serious physics arguments) 4>{t) = 1/ [1 + tan (nt/2)} 
(t < 1) and E cut = 3 x 10 10 GeV. 

Fig. 2 shows the energy dependence of the Z factors, Z Vix (E,x) (solid curves) 
and Z Vji (E, x) (dashed curves) for various depths, from x = to x = x e (where 
x^ ~ 1.1 x 10 10 g/cm 2 is the depth of the earth along the diameter), for the 
initial spectra (9) calculated with 7 = 0.5, a = l(a),7 = l,a = 0.5 (b), 
7 = 1.5, a = 0.5 (c) and 7 = 2, a = 1 (d). In all cases we used E = 1 PeV. The 
calculations were made in the fourth order of the iteration procedure described 
in sect. 2. For all the spectra under discussion, for 10 GeV < E < 10 10 GeV 
and < x < x®, the maximum difference between Z^(E, x) and Z^(E, x) is 
about 4%; the value Z^/Z^ - 1 is less than 2 x 10~ 3 and Z^/Z^ - 1 is 
less than the precision of the numerical integration and interpolation (about 
10~ 5 ) adopted in our calculations. After the tests with many models for the 
initial spectrum, we conclude that the convergence of the algorithm is very 
good and that even the first approximation, Z^(E,x), has an accuracy quite 
sufficient for the majority of applications of the theory. 

As it is clear from fig. 2, the shape of the Z factors is very dependent from 
the initial spectrum. This is a positive fact for neutrino astronomy, since it 
gives, at least in principle, the possibility to reconstruct the initial neutrino 
spectrum from the measured energy spectrum and angular distribution of the 
neutrino induced muon events in a neutrino telescope. 

At comparatively low energies (except for unrealistically hard spectra like 
the one used in fig. 2. a), the Z factors for antineutrinos exceed the ones for 
neutrinos. Considering the inequality X-p^(E) > \ U ^(E), one can conclude that 



A^(E,x) > A^(E,x) 
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E (GeV) E (GeV) 



Fig. 2. Z factors, Z„ (E, x) and Z v ^{E,x) vs energy for the initial spectra (9), 
calculated with four different sets of 7 and a and with Eq = 1 PeV for depths 
x = x^/k [k = 1,2,3,5,10,20,50 from bottom to top] and x = (the largest Z 
factors) . 

for any depth. In the multi-PeV energy region and above, the Z factors (and 
effective absorbtion lengths) are identical for and V^. The difference be- 
tween the shapes of Z v (E,x) and Z Vti [E,x) is almost depth-independent and 
becomes more important for steep initial spectra. This behavior may be under- 
stood from an analysis of the shapes of the total cross sections and regeneration 
functions for and (fig. 1). 

At any fixed energy, the Z factors monotonically decrease with increasing 
depth and the inequality Z V (E, x) < Z®(E) takes place for any x > 0. This ef- 
fect leads to significant decrease of the neutrino event rates in comparison with 
those estimated in the "standard" approximation Z u ps the latter works 
only at low energies, when the shadow effect is by itself small (that is when the 
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medium is almost transparent for neutrinos). Although these conclusions were 
derived from particular models for the initial neutrino spectrum, cross sections, 
and medium, actually they are highly general and model-independent. 

In fig. 3 we present the penetration coefficient, exp [—x/A v (E, x)], in the earth 
for muon neutrinos with initial spectrum (9) calculated with 7 = 0.7 and a = 
( "quasi-power-law" spectrum). The results are presented as a function of E 
for several nadir angles ($) in fig. 3. a and as a function of d for several values 
of E in fig. 3.b. 




Log(E) (E in GeV) nadir angle in degrees 



Fig. 3. Neutrino penetration coefficient in the earth for the quasi-power-law initial 
spectrum with 7 = 0.7 as a function of E at fixed 1? [0° to 90° from bottom to top 
with steps of 10°] (a) and as a function of 1? for E = 10 fe GeV [k = 3, 4, . . . , 7 from 
top to bottom] (b). 

To evaluate the depth, function of 1?, we used the density profile of 

the earth, p(L), given in ref. [13]. The kinks in fig. 3.b are due to the layered 
structure of the earth. 
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A Monochromatic initial spectrum 



Let us consider an initial "spectrum" of the form S(E — E ) with a fixed pa- 
rameter Eq. In this appendix, we will show how this monochromatic spectrum 
transforms at a depth x in a medium. Let us denote the transformed spectrum 
by G U (E ; E,x). This function must satisfy eq. (1) and simple considerations 
suggest the following ansatz: 



G u (E ;E,x) 



5{E - E ) + 9{E ° E E) ME ; E, x) 



-x/\ v (E ) 



, (A.l) 



where the term with 5 function describes absorption of the initial ( "parent" ) 
neutrinos of energy Eq and the next term - the creation and propagation of 
secondary neutrinos with energy E < E . Substituting eq. (A.l) into eq. (1) 
yields 



dM E o;E,x) 1 



dx 



2/0 



$ v (y, E)^ u (E ; E y , x) dy + Q, v (E, E ) 



+V u (E,E )^ u (E ;E,x), ip v (E o ;E,0) = 0, (A.2) 



where Q u (E, E ) = (1 - y ) (y , E), y = 1 — E/E < 1, and V V (E, E ) is 
defined by eq. (7) 

Let us seek the solution to eq. (A.2) in the form 



X X 

Vv(£ ; E, x) = Q u (E, E ) J exp J 



dx" 



1 



L„(E ;E,x) X U (E ) 



l f C u (E ;E,x") 

1 - Z V (E Q ; E,x) 
Xu{E) 



dx' 



(A.3) 



(A.4) 



with Z u (E ;E,x) an unknown positive function. After direct substitution of 
eqs. (A.3) and (A.4) into eq. (A.2) we have 

yo 

Z U (E ; E, x)^(E ; E, x) = j $ v (y, E)^ u (E ; E y , x) dy. (A.5) 



Clearly, Z u {E ;E,x) -> and vp u {E ;E,x) -> x$ u {0, E )/X U {E ) as E -> E 
for any x. 
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The new "Z factor", Z v {Eq\ E , x) , can be found from eqs. (A.3-A.5) by an 
iteration algorithm similar to the algorithm described in sect. 2. Putting Z v = 
as a zero approximation we have 



Vt v (E, E Q ) 
X V (E)V V (E,E ) 



e xV v (E,E ) _^ 



(A.6) 



and subsequently 

yo 

Z£\E Q -E,x) = J Q v {y,E) 



(E ;E y ,x) 
4 0) (E ;E,x) J 



dy. 



(A.7) 



The next steps of the algorithm are quite obvious so there is no need to write 
out the corresponding cumbersome formulas here. 

Let us briefly sketch the leading approximations for ip u and Z v , since they 
contain the main features of the exact solution. As it is seen from eq. (A.6), 
for every E < Eq there is a depth 

X ^E , E) = - -i — In ^4^1 



at which the flux of neutrinos of energy E reaches the maximum. The function 
x*(E ,E) increases when E decreases and tends to the minimum, X U (E ), as 
E — > Eq. At any finite depth, secondary neutrinos "remember" about their 
parents (the E Q dependence does not disappear with increasing depth). Due to 
the nontrivial shape of the regeneration function $^ (see fig. l.b), the spectrum 
of secondary neutrinos is rather complex and transforms fast with depth. 

For x <C A(£'o), the function ip^ behaves as xQ u (E, E ) /\ V (E). Therefore 

yo 

o 



Vt v (E y , Eq) \ V (E) 
Vt v (E, E Q ) X v {Ey) 



dy. 



Taking into account that X V (E) > X U (E ) for E < E (see footnote 4), we get 
the asymptotic behavior of Z^ for x — > oo: 



yo 



z£\Eo;E,x)~ Q v {y,E) 



Q v (E y , Eq) 

&v (E, Eq) 
u 

x exp [—xV u (E, E y )} dy 



\ V {E) - X u {Eq) 



X U {Ey) 

0. 



X„(E )_ 
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With the function ip u (Eo; E, x) in hand, we can obtain the solution to the 
transport equation (1) for any initial spectrum F®(E). Indeed, multiplying 
eq. (A.l) by F®(E ) and integrating over E , we have 



F v (E,x) = J F°(E )G u (E ;E,x) dE 



OO 



d/m*W dE ° . (A.8) 
E 



The first term on the right side of eq. (A.8) describes neutrino absorption 
and the second the neutrino regeneration due to energy loss through the re- 
actions vT — > vX. Eq. (A.8) is in fact equivalent to eq. (2) but, when the 
function ip u (E Q ; E , x) is known, eq. (A.8) becomes much more convenient for 
calculations because ipv is independent from the initial spectrumF]. Due to the 
mentioned equivalence, we can get a useful representation for the Z factor in 
terms of the function i/j u : 



Z v {E,x) 



Xu{E) 



In 



x 



r) v (y,E)ip v (E y ;E,x)e 



-xV v {E,E y ) 



dy 



(A.9) 



It should be noted that the Z factor calculated in the n-th approximation 
using the algorithm (5-6) agrees only numerically rather than analytically 
with that calculated from eq. (A.9), using the iteration algorithm for ip u . In 
particular, substituting ij) v = ip^ into eq. (A.9) yields 

= Z®(E,x), 



Z u (E,x) = ^^-\n 



x 



1 + 



xZ^(E,x) 
K(E) 



where ZjP(E,x) is defined by eq. (8). Thus 



Zi l \E,x) = Zl 1 \E, 



x 



2XJE) 



+ . . 



< 1. 



However, the Z$P(E,x) can be approximated by Z^\E,x) with very good 
accuracy, because xZ^ (E, x) / X V (E) <C 1 in most cases of interest for neutrino 
astrophysics. 

6 However, eq. (A.8) has one evident technical drawback. To use it, one must cal- 
culate 3-dimensional arrays which are hard to interpolate due to the very strong 
dependence of ip v and Z v from their arguments. From this point of view, the algo- 
rithm described in sect. 2 is of course simpler. 



14 



B Neutrino transport equation with a source function 



Here, we briefly show how to take into account the contributions from pro- 
duction of neutrinos through reactions v[F — > u^X (£ ^ £') or the reaction 
chains mentioned in the introduction in the case when these may be treated 
as corrections to the principal solution described in sect. 2 and appendix A. 
Clearly, the problem reduces to the transport equation (1) with a source func- 
tion S U (E, x) on the right side. In line with our general approach, we will seek 
the solution to this equation in the following formQ 



Jb x 

F s (E,x)= J S u (E,x')exp -J 



l-Z v (E,x") 
Xu(E) 



dx" 



dx' 



(B.l) 



with Z u {E,x) a positive function satisfying the equation 



± 

Zu{E,x) = J r] S (y,E;x)$ v (y,E)dy, 



(B.2) 



where we introduced 



Vs{y,E;x) 



F s {E y ,x) 



F s (E,x)(l- y y 



It is easy to verify that Fs(E, x) ~ xS u (E, 0) as x — > 0. Therefore, 

0) 



Vs(y,E; 0) 



s u (E,o)(i- y y 



and this function is assumed to be finite for any E and y. 

The algorithm for the solution to eqs. (B.l), (B.2) is quite obvious: putting 
Zf) = yields 



x 

FP(E, x) = J S U {E, x - x') e^'/M^) dx ' ; 
o 

i 

Zjp(E,x) = J V { °\y,E;x)<S> u (y,E) dy, 



7 We suppose Fs(E,0) = as the boundary condition. 
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etc. The formal question about the finiteness of the involved integrals over y 
is closely related to the very difficult problem of the asymptotic behavior for 
the uN inclusive and total cross sections as E — > oo. This problem is beyond 
the scope of this study, but we can avoid it introducing a cutoff y = 1 — E/E 
(with E 3> E) as the upper limit of the integrals. The reason for such a cutoff 
is in the fact that any physical source function, S u (E,x), must exponentially 
vanish as E — > oo (cf. footnote 3). 



C Neutrino-nucleon cross sections 



According to ref. [13], the inclusive differential cross sections for the reactions 
v^N — > pTX (CC) and v^N — > v^X (NC), where N is an isoscalar nucleon, in 
the renormalization-group-improved parton model are of the form (E ^> M^) 



da C c(y, E) _ 2G 2 F M N E } A(x, Q 2 ) + (1 - y) 2 B(x, Q 2 ) 



dy vr J (Q 2 /M 2 W + 1) 2 

l 

.^.^ ^ dx. 



da NC (y,E) _ G 2 F M N E r A (x,Q 2 ) + (l-y) 2 B (x,Q 2 ) 
dy 2tt J (Q 2 /M 2 + l) 2 



Here Gf is the Fermi constant, M N and M w (M z ) are the nucleon and W- 
(Z-) boson masses, respectively, Q 2 = 2M^xyE, is the squared invariant mo- 
mentum transfer between the incident and outgoing lepton, and x is the usual 
Bjorken scaling variable. For the V^N cross sections, the structure functions 
A, B, A and B in the above formulas should be substituted for the functions 
A, B, A and B , respectively. All these are 



A=^(u v + d v + u s + d s ) + s s + b s , 

A = \ ( u s + ds) + s s + b s , 
B=^(u v + d v + u s + d s ) + c s + t s , 
B = ^(u s + d s ) +c s + t s , 



A o = \ ( L l + L l) K + d v + u s + d s ) + - (R 2 U + R 2 ) (u s + d 9 ) + C, 
A = - (L 2 a + L 2 ,) (u, + d s ) + - (R 2 U + R 2 ,) (u v + d v + u s + d a ) + C, 
Bo = t (R 2 u + R 2 d ) (u s + d s ) + - (p u + L 2 ) (u v + d v + u s + d s ) + C, 
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Bo = 2 (R 2 u + R 2 d) K + d v + u s + d s ) + - (hi + L 2 d ) (u s + d s ) + C, 
where 

C = (L\ + Rl) (s s + b s ) + [l\ + i$) (c s + t s ) , 

= |sin 2 #iy, R u = —2Rd, La = Rd — 1 and L u = R u + 1 are the chiral 
couplings, and 6w is the Weinberg angle. The functions u VyS = u VtS (x,Q 2 ), 
d v , s = d v ^ s (x,Q 2 ), etc. are the distributions of corresponding valence (v) and 
sea (s) quark flavors in a proton. For all constants involved we used the stan- 
dard values from ref. [16]. 
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